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Abstract. Higher-form flux that extends in all 3+1 dimensions of spacetime is a 
source of positive vacuum energy that can drive meta-stable eternal inflation. If the 
flux also threads compact extra dimensions, the spontaneous nucleation of a bubble of 
brane charged under the flux can trigger a classical cascade that steadily unwinds many 
units of flux, gradually decreasing the vacuum energy while inflating the bubble, until 
the cascade ends in the self-annihilation of the brane into radiation. With an initial 
number of flux quanta Qq > N, this can result in N efolds of inflationary expansion 
while producing a scale-invariant spectrum of adiabatic density perturbations with 
amplitude and tilt consistent with observation. The power spectrum has an oscillatory 
component that does not decay away during inflation, relatively large tensor power, 
and interesting non-Gaussianities. Unwinding inflation fits naturally into the string 
landscape, and our preliminary conclusion is that consistency with observation can be 
attained without fine-tuning the string parameters. The initial conditions necessary 
for the unwinding phase are produced automatically by bubble formation, so long as 
the critical radius of the bubble is smaller than at least one of the compact dimensions 
threaded by flux. 
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1 Introduction 

Many theories, including string theory [1], predict the existence of metastable phases 
with positive vacuum energy. Such phases drive eternal inflation and are expected to 
exponentially dominate the global volume of the universe. If so, our Hubble volume 
must be contained inside a pocket or bubble embedded in the parent false vacuum. 
To be consistent with observational constraints on spatial curvature, the radius of our 
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Dramatis Personas 


D 


Total dimension of space-time 


F,p 


F is a p + 2-form field strength under which a j9-brane is charged 


Q(Qo) 


Number of F flux quanta (prior to bubble formation) 


M 


A compact (D — 4)-manifold 


H,N 


Hubble constant of the 4D space-time, number of efolds of inflationary expansion 


z bi z 


Radius of the brane in the compact dimensions, effective 4D inflaton field 




Velocity v — z, Lorentz factor, rapidity, and tension of the brane 




Speed of sound for perturbations 6z, cf. Sec. 3.2 


I 


Circumference of compact directions threaded by F in which the brane expands 


d 


Circumference of compact dimensions transverse to F or wrapped by the brane 


c 


FRW curvature perturbation, cf. Sec. 2.1 


V ( (k) 


Dimensionless power spectrum of ( (3.5) 


n(p s ) 


Number (energy) density of strings produced per brane/anti-brane collision 


X 


Time averaged friction due to string production, (3.14). 


f 


Time averaged pressure due to string production, see (3.11) 



bubble must be at least 10 times larger than the present-day Hubble length, meaning 
it must have undergone N > 60 efolds of inflationary expansion after it formed. 

In this paper we present the details of "unwinding inflation" [2] , a novel mechanism 
by which an eternally inflating metastable false vacuum can transition via charged 
brane bubble formation to flux discharge cascade [3] that mimics slow-roll inflation. 
The vacuum energy is initially only slightly reduced by the formation of the bubble, but 
then steadily "unwinds" over time. Inflation ends and reheating occurs with the self- 
annihilation of the brane into radiation once most or all vacuum energy is discharged. 
Without fine-tuning of either the parameters or the initial conditions, this phase can 
drive 60 or more efolds of expansion, solving the curvature problem and generating a 
scale-invariant spectrum of perturbations. 

Unwinding inflation occurs when the vacuum energy of the parent phase is at 
least partially supplied by electric-type flux that extends in the 3+1 dimensions of 
spacetime, and in addition wraps at least one compact extra dimension (hence, a five 
or higher form flux). Under these circumstances, a brane bubble that reduces the 
flux by one unit of the brane's charge can nucleate [4]. If the bubble is localized on 
the compact space, it will expand, rapidly reaching relativistic velocities, wrap around 
the compact dimension, collide with itself and initiate a flux cascade that repeatedly 
discharges the flux one unit at a time [3] . The effective four-dimensional inflaton scalar 
z — the radius of the bubble in the compact dimensions — increases steadily with time, 
as in conventional slow-roll models. If the size of the extra dimension is I and the 
Hubble constant is H, roughly (Hl)^ 1 units of flux will be discharged per Hubble time. 
We believe compactifications with HI > 1 are difficult to achieve, and therefore to 
attain N efolds of inflation probably requires Qq > N, where Qo is the number of 
initial flux units. 

A heuristic analogy is that the flux is a rubber sheet that wraps multiple times 
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around some compact cycles, as well as extending in the 3+1 large dimensions of 
spacetime. The nucleation of the brane bubble is the spontaneous appearance of a 
spherical hole (bounded by brane) in one layer of the sheet. The total energy of the 
brane when it appears is equal to the total energy cut out of the sheet. Once it appears, 
the tension of the sheet pulls on the hole and causes it to expand in all the directions 
the sheet extends in. As the hole expands around the compact cycle and overlaps itself, 
it unwinds more and more layers of the sheet. 

This mechanism uses ingredients (branes and flux) found in all string compact- 
ifications, and as such fits naturally into the string theory landscape [1]. The cos- 
mological constant problem — the requirement that inflation ends at or close to zero 
vacuum energy — is solved in theories with sufficiently large numbers of eternally inflat- 
ing phases [5]. String or M-theory compactified on a manifold with length scale O(10) 
times larger than the fundamental length allows for a more than sufficient number [6], 
and also allows Qo ~ 100 as required for unwinding inflation. The two hierarchically 
different scales involved in bubble formation followed by slow roll that appear unnat- 
ural from the low-energy point of view [7] (the thinness of the bubble wall versus the 
flatness of the inflationary plateau) arise naturally: the bubble wall is a D-brane and 
therefore very thin, while the flatness of the inflationary plateau follows if Qo > N. 

Unwinding inflation has the virtue that it sets up its own initial conditions — 
eternally inflating false vacuum states are exceptionally powerful attractors [8], and 
unwinding inflation initiates spontaneously from them. At the same time, this model 
reheats homogeneously and isotropically and explains how the bubble can become 14 
Gyrs in size, thereby solving the problems that derailed Alan Guth's original model 
of "old" inflation [9]. Unlike in old inflation, the rate of bubble nucleation can be 
arbitrarily small — a bubble will nucleate eventually, and only one is needed. 

1.1 Observational consequences 

Several observational tests of the eternally inflating multiverse have been proposed, 
among them searching for the effects of cosmic bubble collisions (see [10] for a recent 
review) and measuring non-zero spatial curvature [11, 12]. The weak point in these 
tests is the sensitivity to the duration of slow-roll inflation after the bubble forms. Long 
inflation in the bubble erases all signatures of such "initial state of the bubble" effects 
with exponential efficiency. While there are reasons to believe inflation was short [13], 
that conclusion is weakened by the lack of understanding of the origin of the slow-roll 
phase. 

By contrast, unwinding inflation predicts a set of characteristic features that do 
not inflate away. Fundamentally, this is because the compact dimensions remain stable 
and small, but nevertheless play a continuous and key role during inflation. Brane col- 
lisions occur periodically as the bubble expands around the compact dimensions and 
intersects itself, producing open and possibly closed strings. This string production 
contributes an oscillatory component to the power spectrum of perturbations. The in- 
flaton potential arises primarily from the background flux, but it too has an oscillatory 
component due to brane-anti-brane interactions in the compact directions. 
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The presence of Lorentz invariant (Dirac-Born-Infeld) kinetic terms for the brane 
and the possible presence of additional light scalars (describing the transverse position 
of the brane in the compact dimensions) lead to non-Gaussianity. Unlike most pre- 
vious attempts to realize inflation in string theory, the scale of inflation is high, and 
unwinding inflation naturally predicts tensor modes with an amplitude that can be 
observed in the near future. 

The exact details of our observational predictions are necessarily tentative, be- 
cause we do not yet have a realization in a fully stabilized string compactification, and 
by the same token do not know the geometry of the compactified manifold. We will 
attempt to elucidate to what extent our predictions are generic and independent of the 
details of the compactification, and where they may break down. At least in principle, 
confirmation of the predictions of this model could provide evidence for the eternally 
inflating multiverse and the presence of cosmic bubble collisions (albeit in a compact 
dimension), probe the geometry of the extra dimensions, and provide observational 
support for string theory. 

1.2 Relation to previous work 

A short, self-contained description of this model can be found in [2]. Discharge of 
higher- form flux by branes was first considered in [4] , and flux discharge cascades in [3] . 
The idea of using a bubble collision in a compact dimension to reheat homogeneously 
and isotropically was first proposed in [14], and the possibility of scalar cascades was 
mentioned in [15, 16]. Fluctuations of bubble walls were considered in [17-19]. Another 
model of inflation that uses compact extra dimensions to extend the field range is 
[20], related 4D effective field theories were studied in [21-23] and the idea of using 
a relativistic brane for inflation with a Dirac-Born-Infeld (DBI) action was proposed 
in [24, 25]. The effects of particle production on the spectrum of perturbations was 
considered in [26-28]. Various models have been considered in the past that utilize 
the attractive potential between a brane and anti-brane for inflation (e.g. [29]). Such 
models typically need warping to make the potential flat enough (e.g. [30]). Unwinding 
inflation does not, because the inflaton potential arises instead from background flux 
(brane/anti-brane interactions are a small, periodic perturbation in the potential). It 
also evades the exponential rarity of post-tunneling slow-roll described in [31], because 
the tunneling creates the conditions necessary for slow roll. A model using multiple 
quantum tunnelings in place of slow roll (in contrast to unwinding inflation, which is 
classical after an initial bubble nucleation) is [32] . After this paper was completed and 
[2] appeared, [33] was posted, which overlaps with [2, 3], and this work. 

2 Background and basic mechanism 

Higher-form fluxes that extend in the 3+1 large dimensions of spacetime and in any 
number of compact extra dimensions contribute to the 4D components of the stress 
tensor like vacuum energy. In string theory (or any other theory with several extra 
dimensions) the different ways of threading flux through compact cycles gives rise to 
a large landscape of metastable vacua. With the mild conditions that the compact 



-4- 



dSi compact manifold Ai 
, * , , A , 

0123 45 6 7 8 ••• 

-k -k -k -k -k -k 

V V 

Flux fills these directions Transverse dimensions 

Figure 1. The setup for unwinding inflation in the case p = 4. 



dimensions be larger than the string length by an 0(1) factor and that there are 
(9(100) distinct compact cycles, the number of possible vacua can greatly exceed 10 120 , 
thus providing the necessary conditions for an anthropic solution to the cosmological 
constant problem [5, 6]. 

Typical vacua of this type have a large vacuum energy and innate very rapidly, and 
are therefore expected to dominate the overall volume of the universe. Bubbles of other 
phases form inside these rapidly inflating regions via a variety of phase transitions such 
as the tunneling of metastabilized scalar moduli from a minimum and the discharge 
of a unit of flux by the nucleation of a spherical bubble of brane [4] . Both transitions 
produce bubbles containing FRW cosmologies, and additional inflation is necessary to 
solve the curvature problem. 

Our mechanism occurs when a bubble of brane (that initially discharges one brane 
charge unit of flux) appears and is smaller than any of the compact directions it extends 
in (namely, the directions the flux threads). Under these conditions there is a "flux 
discharge cascade" [3], where a (p+2)-form electric flux threading at least one compact 
dimension can "unwind" , repeatedly discharging in a cascade triggered by the quantum 
nucleation of a bubble of charged brane, and hence steadily decreasing the effective 
four dimensional vacuum energy. 

In the case of a top (D) form flux, the brane is co-dimension 1 (a domain wall) 
in the space. In that case the flux everywhere inside and outside the brane bubble is 
constant, but the interior flux is reduced relative to the exterior by one unit of the 
brane's charge. The non-zero flux at the bubble wall exerts a force on it, causing 
it to expand in all directions after it appears. The bubble expands freely in the 3+1 
dimensions of spacetime, but in the compact dimension (s) it wraps around and collides 
with itself. Ignoring self-interactions for a moment, it passes through itself, and in the 
overlap region discharges the flux by two units. This region expands and, after a second 
wrap, forms a region with three units discharged, etc. (Fig. 2). 

The necessary ingredients for unwinding inflation are (Fig. 1): 

• A (p + 2)-form field strength F with p > 3, and p-branes that are electrically 
charged under F. 

• A D = 4 + q dimensional spacetime dS± x M. q) l where dS& is 4D de Sitter 
spacetime and Ai q is a stabilized compact g-manifold with q > p — 2 > 1. 

• Qo ^> 1 units of F flux threading dS^ and a p — 2 cycle in Ai q , supplying the 
dS± vacuum energy. 

1 One could generalize this and consider a warped product, but we will not do so here. 
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Figure 2. The mechanism of flux discharge cascade on dS^ x Si, where z ~ z + I is the 
coordinate on the Si- The amount of flux is indicated by Q, Q — 1, . . ., while the dashed arrows 
represent the direction of the electric force and the velocity of the branes. The squiggly lines 
indicate strings stretched between sections of brane; their mass depends on the separation 
and changes with time. The figure is not to scale; usually I < 1/H,R where R is the radius 
of curvature of the branes. 

An observer located at a point in the spacetime would encounter a series of brane 
walls that sweep across her location at regular intervals, and a flux that is constant ex- 
cept when a wall crosses her position, after which it decreases by one unit. An observer 
unable to resolve distance or time scales of order the size of Ai would simply observe 
a steadily decreasing flux. Because the flux contributes positively to the effective vac- 
uum energy of the 4D spacetime, during the cascade there is a gradual decrease in the 
Hubble constant of the de Sitter — just as in ordinary slow-roll inflation. After part or 
all of the flux is discharged, the remaining effective 4D vacuum energy can be positive, 
negative, or zero, depending on the stabilization mechanism and any additional fluxes 
or vacuum energy. For our purposes, we will assume the 4D vacuum energy at the end 
of the cascade is close to zero, or at the value that will result in nearly zero vacuum 
energy after GUT or standard model phase transitions that occur later in the evolution 
of the universe. 

The reduction in flux during the cascade can also lead to a change in the geometry 
of Ai, typically reducing its overall volume. However if the flux is not the primary 
element that stabilizes Ai this change is small, and we will neglect it (we will comment 
more on this in Sec. 4). 

2.1 Prototype model 

The simplest version of unwinding inflation has D = 5 and p = 3. In this version 
the initial spacetime is dS± x Si, where the vacuum energy of the dS± is supplied in 
part by Qo units of initial F$ flux. Stabilizing the Si requires an additional ingredient, 
for example the Casimir energy of several bosonic and fermionic fields [34]. In Sec. 4 
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we will describe stabilization in more detail, but for now will simply assume the Si is 
stable. The objects charged under F 5 are p = 3-branes, with charge e and tension a. 
With Q units of flux, one has F§ = Qe. In some cases, higher dimensional versions of 
unwinding inflation can be dimensionally reduced to this 5D version. 

Instanton: The nucleation of brane bubbles that discharge the flux is governed by a 
solution to the Euclidean signature equations of motion [4]. The Euclidean signature 
version of the spacetime has metric 

ds% = H- 2 (d£ 2 + sin 2 £dttl) + dz 2 , (2.1) 

where z ~ z + I is the coordinate on S\ and d£l\ is a 3-sphere. 

Typically, the dominant instanton for decay of a false vacuum state has the max- 
imal symmetry possible. We will assume the brane is thin, and that the instanton 
depends only on £ and z in accord with the symmetry of the initial state (so that the 
bubble is spherical in the dS± directions). We are interested in the case where the 
maximum size Az of the instanton in the z direction satisfies Az < I, so the periodic 
boundary conditions do not affect the solution (at least in the thin- wall limit). Finally, 
when the initial number of flux units Qo S> 1 the gravitational backreaction of a single 
bubble is small and can be ignored. 

With these assumptions the instanton solution is fully characterized by the lo- 
cation of the wall z = ±Zi,(£), where ± refers to two symmetric halves and we have 
chosen z = £ = as the center of the bubble. To find z b (^), one should minimize the 
action Se = —k f v ^/Qe + ct f dv ^/gE {induced), where k is the difference in energy density 
on the two sides of the wall, a is the tension of the wall, and V is the volume enclosed 
by the bubble. 

Parametrizing the position of the wall by z = ±2&(£), the Euclidean action is 
S E =H~ A I dz I dlV^sin 3 ^ - kG(z + z b )G(-z + z b ) 



1/2 

1,1- \ - 

r2 



dz x 



+ 2a5(z -z b ) l + H y sji 
-An 2 H- 4 J sin 3 £ I -Kz b (() + a 



1 + 



(2.2) 

1/2" 



Extremizing this action results in equations of motion that can be solved ana- 
lytically for dz b /d^, with the integration constants fixed by the requirements of finite 
action and smoothness: 

dz b z(8 + 9cos£ — cos3£) (2 3) 

d£ ~ Hy/-(8 + 9 cos £ - cos 3£) 2 + lU(aH/n) 2 sin 6 £ 

This solution describes a bubble with spherical topology in five dimensions. The shape 
is oblate; the z direction is distinguished from the Euclidean de Sitter directions because 
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it is flat, and the coordinate extent in z is somewhat less than it would be in 5D 
Euclidean space. 

We can determine the Lorentzian signature evolution of the bubble after it nucle- 
ates through the analytic continuation £ — > iHt, dQ 3 — > idH^. In flat space, the bubble 
would expand and accelerate indefinitely with constant proper acceleration, reaching 
a gamma factor 7 = z/2R after expanding by a distance z [3]. Instead, because of 
Hubble friction the wall approaches an asymptotic velocity v < 1 in the ^-direction 
[14] 

]im^ = v = 1 1 (2.4) 

M y/l + {ZoH/k) 2 sj I + (3R H/4) 2 

where Rq = 4a /k would be the radius of the bubble in five dimensional flat space 
(H = 0). Additional sources of friction (for instance from string production due to 
collisions) may reduce v even further. 

Because their charges are opposite, it is useful to think of the half of the bubble 
at z = +Zb as brane, and z = —Zb as anti-brane. Unwinding inflation eventually ends 
when the brane annihilates with an image anti-brane. For various reasons detailed 
below, the brane position Zf, will become inhomogeneous with perturbations Szb(x,t). 
Assuming that the unperturbed brane will annihilate with its image after Q t units of 
flux are discharged, reheating occurs when 

z b (t) + Sz b (x, t) = vt + 5z b = Qtl/2. 

Solving for t yields 

t = Q t l/2v - 5z b /v = t + St. 

The curvature of this hypersurface is a(t + St) = a(t ) + a(t )St, and so the curvature 
perturbation is 

( = Sa/a = HSt = HSzb/v = HSz b /zb, 
where H and £& are evaluated at horizon crossing as usual. 

Effective action: After formation of the bubble, the brane/flux action in Lorentzian 
signature is 



S = H- 3 J dz J dH 3 dt smh\Ht) (-2aS(z - z b )y/l - (dz b ) 2 - 



(2.5) 



where dH 3 = sinh 2 p dp (KI2 is the measure on a unit hyperboloid. Gauss' law [4] 
requires that the flux changes across the brane by one unit of the brane charge: 



F i _ ..5 



5! 



irQ 2 = /1 5 Q + Yl ~ z » + ^ ~ & ( z + Zb + ^ 



Here p 5 ^ 2 is the charge of the brane, Qo is the number of flux units prior to the bubble 
nucleation, and the sum arises due to the periodicity of z. From this, one can see that 
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the jump in vacuum energy across a brane separating Q from Q — 1 units of flux is 
K(Q) = / iS(Q- 1/2). 

As usual for cosmic bubbles that nucleate due to first-order phase transitions, 
the bubble's walls and interior are naturally described using the negatively curved 
spatial slicing of (2.5). In these coordinates, the energy density inside the bubble is 
homogeneous and isotropic. The spatial curvature of the bubble universe is = 
(aH)~ 2 « e~ 2N after N efolds of inflation, and today would be flk ~ e ~ 2 ( N - N *)^ where 
N* ~ 60 and N is the total number of efolds of inflation after the bubble forms. Since 
observation constrains Qk ^ -01 [35], unwinding inflation requires N > iV* + 3. 

After a few efolds of expansion, the radius of curvature of the bubble R(t) ~ 
Roe Ht stretches to super horizon scales and Qk is exponentially small. From then on 
during inflation, to a good approximation the bubble can be treated as a flat, parallel 
brane/anti-brane pair separated by a distance 2zb in the compact dimension, and the 
open slicing used in (2.5) can be replaced with the flat slicing (see Fig. 2). We will 
neglect corrections due to spatial curvature for the remainder of this paper. 

The action (2.5) can be integrated over z to obtain a 4D effective action: 

S = J dtd 3 xe 3Ht (-2a^l-(d t z b y + e-™t{dtz b y - V(z b fj , (2.6) 

where we have approximated the open slicing in (2.5) with flat slices, and V(zb) is a 
piecewise-linear interpolation of a quadratic that is most easily expressed in terms of 
its derivative (see Fig. 3): 

Z - - (<* - \- 

where [...] denotes integer part. 

This expression has a very simple physical interpretation: —V'(zb) is the pressure 
on the brane, which is proportional to the (integer) number of flux units at its location. 
Every time [2zi,/l] increases by one, an additional collision has occurred and the flux 
(and hence the pressure V) is reduced by one additional unit of the brane charge /x 5 . 

3 Four dimensional effective description 

In the previous section we derived a 4D effective action (2.6) describing the simplest 
version of unwinding inflation. As we will see, other compactification geometries pro- 
duce 4D actions that differ in some details, but all of them share certain features 
in common. In this section we will derive the cosmological predictions of unwinding 
inflation in terms of a few parameters of the 4D action. 

The effective 4D dynamics of unwinding inflation can be described by the following 
ingredients. 
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Figure 3. Left panel: The pressure V from (2.7), constant between collisions and changing 
in steps at each collision. Right panel: The potential V around the zero flux minimum; the 
smooth line is the quadratic approximation. 

• The inflaton scalar field Zb (corresponding to the radius of the brane bubble in 
the compact dimensions that are threaded by flux). 2 

• D — p — 2 additional scalars b (corresponding to the position of the brane in the 
compact directions perpendicular to the flux). 

• Dirac-Born-Infeld kinetic terms for all the scalars, normalized by the tension of 
the brane a and possibly functions of the fields. 

• A potential V(zb) (arising from Q units of electric flux that force the bubble to 
expand, as well as from the interactions of the charge of the brane with its images 
in the compact directions). 

• A collection of extra degrees of freedom that are produced periodically in bursts 
when the branes collide. 

The bubble nucleation sets up initial conditions in which all fields are homogeneous 
and isotropic on negatively curved (open universe) slices, and V(z b ) ~ Qq is relatively 
large All cosmological consequences of this model can be determined from the 4D 
description. At the level of detail we will present here, the relevant quantities are a, 
V, the average density of produced strings/particles/gravitons p s , and the number of 
scalars along with information regarding how they influence the time of reheating and 
modulate the kinetic terms. 

In general, the reduction to four dimensions produces a Kaluza-Klein-type tower 
of modes with masses m ~ 0(1/ I), arising from excitations of the brane in the compact 
directions it extends in. Since we consider only compactifications with HI < 1, these 

2 It is important to note that Zb is not a fundamental scalar — it is a collective coordinate that does 
not exist prior to bubble nucleation. Its potential V(zb) does not in general have false minima or a 
tunneling transition. 
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modes are heavier than the Hubble scale and do not affect the dynamics of inflation 
significantly. The 4D effective action for z (for clarity from here forward we drop the 
subscript on Zb) is: 



S 




(c.f. (2.6), (A.5)). 



The potential V has contributions both from the unwinding flux and from other 
sources, such as the higher dimensional cosmo logical constant, other fluxes, Casimir 
energies, etc. Assuming the extra dimensions remain stable throughout the cascade, 
these latter contributions are constant. By contrast, the energy density due to the 
unwinding flux scales as Q 2 and therefore depends on z and decreases during inflation. 
As we mentioned above, the 4D vacuum energy must be tuned so as to solve the 
cosmological constant problem. This means the constant term is a small contribution 
relative to the ~ Q 2 piece, and we will neglect it. 

During unwinding inflation when Q > 1, the kinetic energy 2ah / y arising from 
the first term in (3.1) is small compared to V. This must be the case at least early on, 
because the brane nucleation conserves energy, and therefore the energy in the brane 
tension equals the energy released by discharging only one unit of flux. Furthermore, 
Hubble friction limits the Lorentz factor 7, and so it is only when most or all of the flux 
has discharged that the kinetic energy becomes important. Therefore, during much 
of the cascade the 4D Einstein equations are H 2 ~ V, as usual in models of slow-roll 
inflation. 

A probe brane moving on a warped product manifold (where the radius of the dS$ 
depends on position in the compact dimensions) gives rise to an action similar to (3.1) 
(see e.g. [24]), but we are not considering warping here. Instead, in uwinding inflation 
the function h(z) arises because z is the radius of the bubble expanding in a (possibly 
curved) compact manifold. We give some examples of compactifications that produce 
non-trivial h(z) in Appendix A and Appendix B, but for simplicity in the remainder 
of this section we will focus on the case where h(z) = 1: 



The action (3.2) is incomplete in that it only describes the radius of the bubble, 
and - in the string theory context - ignores the other open string degrees of freedom, 
as well as the coupling between open and closed strings (apart from the coupling to 
the background flux, which is taken into account in V). The massless scalar that 
describes the center-of-mass position of the bubble in the compact directions does not 
contribute to curvature perturbations [17], because its value does not affect the time of 
reheating (at least at lowest order). More importantly, strings stretched between the 
brane and its images have a mass that depend on z (see Fig. 2). Because z is changing 
in time these open strings are produced when the brane scatters [36] off its images in 
the compact dimensions, in a process very similar to that of [27]. Furthermore, closed 
strings may be produced by Bremsstrahlung. Here, rather than explicitly including 




(3.2) 
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Figure 4. From left to right: The brane position z, the Lorentz factor 7, the oscillations 
around the smooth approximation of the potential, zoomed in to the region corresponding to 
the CMB quadrupole. All plots use the set of parameters: g s = .01, I = 19.7/m s , d = 2/m s , 
Qo = 400, for a wrapped (p = 4)-brane expanding on an Si, see Sec. 5. As is apparent from 
the plots the oscillations in 7 and V are very small in this case, but the power spectrum may 
have larger oscillations depending on the degree of string production. 



these modes, we will incorporate this effect by including the average amount of particle 
production per collision p s (see [27, 28]). 

3.1 Background evolution 

The potential energy is the energy in the flux: V oc F 2 oc Q 2 , averaged over the 
compact directions. The specific form V(z) depends on the compact geometry, but 
in all cases it decreases with increasing z as Q discharges, and includes an oscillating 
component arising from interactions of the brane with itself as it wraps the compact 
dimensions (c.f. Fig. 3, Fig. 7). In several examples 3 , V(z) ~ Q 2 ~ {Qo — z/l) 2 + small 
oscillations, where I is the typical lengthscale of the compact cycle the flux threads. 
The amplitude of the oscillations is suppressed by at least 1/Q, since it arises from 
interactions of the brane with itself rather than the Q units of background flux. These 
examples are similar to m 2 2 inflation, although with a DBI kinetic term and extra 
degrees of freedom. 

Generically the collision or scattering of the brane with its images is inelastic due 
to the creation of particles or strings. This particle creation will affect the background 
equation of motion (3.11) because it converts some of the kinetic energy of the brane 
into particles, therefore acting as a source of friction. However, in most regimes of 
interest in unwinding inflation this additional friction term is subdominant to Hubble 
friction. 

Ignoring the effects of particle production for now, the background (homogeneous) 
equation of motion from (3.2) is 



3 When the cycle the flux threads is S n for any n, or when it is S n x Af with the size of AT small 
enough that the brane bubble wraps it rather than expanding in it. 
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The combination of the force due to V, the DBI kinetic term, and Hubble friction 
results in z ~ vt with the slow-roll velocity v nearly constant, with 1 — v <C 1 in most 
cases. The small oscillations in V do not strongly affect the background evolution 
(see Fig. 4). 

As a check on our formalism, for the "prototype" model where V = 2k = 2/i 5 (Q — 
1/2), v — z obtained from (3.3) with z = agrees exactly with the result derived from 
the instanton (2.4). 

3.2 Perturbations 

There are two sources of perturbations (both scalar and tensor) in this model. First, 
the finite temperature of de Sitter space leads to tensor and scalar fluctuations in the 
standard way. As we will see, the power spectrum of these perturbations depends 
only on the Hubble constant H during inflation, the effective brane tension a, and 
the velocity z (and h(z) when it is not constant). Second, the collisions of the brane 
bubble with itself as it wraps around the compact dimensions will produce particles 
and/or strings. Determining the rate at which this happens requires a model (see 
Appendix C), but the effect on the power spectrum of perturbations depends at first 
approximation only on the energy density produced as a function of z and its time 
derivatives. Particle production is a random process, and (at least assuming the theory 
is weakly coupled) particle/string production events should be Poisson distributed 
at separations larger than the Compton wavelength of the produced particles. This 
assumption determines the statistics of fluctuations in the rate of particle production, 
and therefore the statistics of the resulting fluctuations in z. 

The brane collisions could also give rise to gravitons, either Bremsstrahlung from 
scattering or via the decay of some other produced particle. It is possible for the 
contribution of tensors produced this way to exceed those produced by standard de 
Sitter fluctuations [37], but we will defer investigating that possibility to later work. 

De Sitter perturbations: We begin by analyzing the model in the absence of any 
string or particle production. With this simplification, the power spectrum for theo- 
ries of the form (3.2) can be computed using standard methods [38]. The late time 
expression for 5z is very simple: 

fc ^^P' (3 - 4) 

As usual, the curvature perturbation is related to perturbations in z by ( = H5z/z. 
Defining the power spectrum by 

(((k)C(k>)) = 2n 2l ^5 3 (k-k'), (3.5) 

we find 
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Incorporating non-trivial h(z) (see (3.1)) multiplies (3.6) by l/h(z), but this is reliable 
only when h(z) and the speed of sound c s vary slowly. In the case of unwinding inflation 
on S2 for example, this adiabaticity assumption does not hold (see Appendix A). 

Because observations determine V( ~ 2 x 10~ 9 , the Hubble parameter H in un- 
winding inflation satisfies H ~ 0.02 a 1 / 4 . Hence if the brane tension a is close to the 
string or Planck scale (as it is in string theory), unwinding inflation is a high-scale 
model of inflation. 

The speed of sound — which is c s = I/7 for (3.1) and (3.2) — cancels out of the 
power spectrum (3.6), although as we will see it contributes to non-Gaussianity. This 
cancellation can be understood most simply as follows. The power spectrum can be 
written as V$ ~ G 2 v V 2 /(c s T), where Gjy is the 4D Newton constant, V is the potential 
energy and T is the kinetic energy. In our case, T ~ 70" = cr/c s , so the speed of sound 
cancels. 

Tensor power: The tensor power spectrum due to de Sitter perturbations with 
action (3.2) is simply [38] 

V h ^-^, (3.7) 

7T 

therefore the tensor-to-scalar ratio is 

\2%TxG N av 2 



H 2 

As we will see in Tab. 1, in models derived from string theory, r is potentially observable 
in the near future. 

Oscillations: The power spectrum (3.6) oscillates due to the oscillations in V ~ H 2 . 
For the "prototype" model (2.6), the amplitude of these oscillations relative to the 
average is very small, ~ Q~ 2 (Fig. 4). For the S 2 (discussed in Appendix A), the 
oscillations in the V are ~ l/Q. 

Tilt: Because v ~ 1 and assuming o is constant, the tilt of the scalar spectrum 
n s — 1 = dlnV /dink arises from the if 4 term in (3.6): 

.H dt H , n , 



During unwinding inflation H ~ y V ~ Q, and Q is constant at least when V'(z) ~ z. 
In this case, 

n,-l« -2/iV* w -0.033, 

where N* = J Hdt = J dH(H/H) —H 2 /(2H) is the number of efolds from the time 
the quadrupole mode crossed the horizon during inflation to the end of inflation. In a 
model with a high reheat temperature iV* ~ 60, and the tilt (3.9) is consistent with 
observation [35]. 

Ordinarily V" contributes to the tilt because it is a mass term for the perturbation, 
and variations in the speed of sound contribute in the form c s /Hc s . In our case these 
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two contributions exactly cancel to first order in slow roll, as is evident from (3.9) and 
can be checked using eq. (35) of [38]. 

A potentially significant correction to (3.9) could arise from changes in a during 
inflation. As we discuss below, if realized in string theory the brane is likely to be 
wrapping a compact cycle, in which case a depends on the volume of that cycle. If 
the volume of the compact cycle changes significantly during inflation, a will change 
as well, therefore tilting the spectrum. However at least in toy models (see Sec. 4), 
the size of the compact cycle changes only slightly as the flux is discharged, making a 
small or negligible contribution to the tilt. 

3.3 Effects of particle or string production 

Each time the brane collides with itself, some fraction of its kinetic energy will be 
converted into particles or strings. The particles produced by any given collision will 
dilute away exponentially on Hubble time scales, but because collisions happen at least 
once per Hubble time, the time-averaged density is not necessarily small. 

Here, we will include the effects of particle production without referring to the 
underlying model. Our treatment is in many ways parallel to [27, 28]. Specifically, the 
field theory model considered in [27] is the non-relativistic limit of the sector of our 
model that describes the coupling of z to the modes of the open strings that are massless 
when the brane and anti-brane coincide. The work of [28] considers dissipation in more 
generality. Our methods differ from these works in several ways, one being that we do 
not assume that dissipative effects are the dominant source of perturbations or friction 
in the background equation of motion. Instead, we will derive a general result that 
applies to both those cases and situations where de Sitter perturbations dominate. 

We will assume that the produced particles are massive and weakly interacting, so 
that the number density n c produced at each collision redshifts as a" 3 (this assumption 
can easily be relaxed), and that they are created instantaneously at times when the 
branes collide. The energy density in the produced particles is 

p s = ^m c ( z , Zi , Zi, ...) n c {zi, z h ...) e- m{t - u) Q(t - U), (3.10) 

i 

where m c is the mass of the string or particle produced by the collision, and 2j, 
etc. are the position, velocity, etc. of the brane at the collision times t*. In string 
theory as well as in the field theory considered in [27], the mass of the produced 
particles/strings grows linearly with z — z iy and we will assume m c oc z — Zi, although 
this can be generalized without difficulty 4 The number density n c depends on the time 
derivatives of z at the collision time, because the rate of change of the mass and the 
kinetic energy at the collision determines the amount of particle production. 

4 The reader might be concerned that this assumption implies that the particles are massless when 
created, but that the modes of stretched open strings do not necessarily have this property, depending 
on their oscillator mode. However, in string theory the string production occurs in a time At ~ 
l/m s <C 1,1/ H, and the typical mass of the stretched strings before they begin to dilute due to 
Hubble expansion is m 2 s /H 3> m s , so corrections to (3.10) from string mode excitations are small. 
See Appendix C for details. 
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We can incorporate the effects of string production and derive the homogeneous 
background equation for z using the continuity equation p = —3H(p + p) as well as 
(3.10) and (3.2): 



= 2j 3 azz + QHaz^ + V'{z) + % + 3Hp s 

at (3.11) 

= 2^az + QHaz-i + V\z) + /, 

where f(z, z, ...) = dp s /dz is the "force" due to string production, and the second line 
follows when m c oc z — Z{. 

The force due to the background flux V 3> / in the early stages of inflation 
(although near the end of inflation when enough flux has been discharged this is no 
longer the case) so the term proportional to / is a small correction to the background 
evolution. Nevertheless we will retain it, as it plays an important role in the dynamics 
of the perturbations. 

Perturbations: The variation in / arises from several sources. First, variations in 
the velocity z affect the amount of string production (we will neglect any dependence 
of / on z and higher derivatives). Second, variations in the time of collision 5ti = Sz/zi 
affect the result. Lastly, string or particle production is a quantum process, and there 
will be random variations Sn in the number density of produced strings. Putting this 
together, the variation in / is 

i i v N ' 



5t l -^n c (3HQ(t-t i )-5(t-t l )) + 5n c -^Q(t-t l )j w Sz + m 2 5n. 

(3.12) 

Here, m\ = dm c /dz, and the overbar and refer to a time average. In what follows, 
we will solve the time-averaged equation, since the equation that results from keeping 
all the terms in (3.12) is difficult to deal with. We will comment on the accuracy of 
this approximation shortly. 

Using (3.12) and (3.11), after time- averaging the perturbed continuity equation 

is 

Sz + 3H(1 + X)5z - e~ 2Ht %5z = , (3.13) 

where we have dropped a term proportional to V", and A is a dimensionless parameter 
controlling the strength of the friction due to the string or particle production in (3.13): 

d f 

A = -^-. (3.14) 

To proceed, we need to evaluate the source term 8n that arises from fluctuations 
in the number density of produced particles or strings. Assuming particle production 
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is a Poisson process in physical space, the power spectrum of perturbations in the 
Fourier transformed number density 5n^ is 

(<%<%,) = ^{2<nf5\k + k') , (3.15) 

where the factor of a~ 3 is needed because k is comoving. Eq. (3.15) implies 

8n % = Vna~ 3/2 X % , (3.16) 

where X^ is a stochastic field satisfying (X^Xp) = (2ii) 3 5 3 :i (k + k'). 

Power spectrum: Equation (3.13) in momentum space describes a damped har- 
monic oscillator with an external source proportional to 5n. The general solution to 
(3.13) is a linear combination of two independent solutions of the homogeneous equa- 
tion plus a particular solution. The latter can be conveniently expressed in terms of 
a Green's function integral. The homogeneous solutions (in terms of conformal time 
r = —e~ Ht /H) are 

<%(r) = C 1>k r v J v {c s kr) + C 2>k r v Y u (c s kr) , (3.17) 

where J u and Y u are Bessel functions, and v = (3 + A)/2. 

At late times (r — > 0) the first term vanishes, so only the second term contributes 
to the power spectrum. To fix the coefficients, we need to specify the initial conditions 
for the mode. For frequencies higher than the energy of the produced particles, the 
particle production should not affect the modes of 8z. Therefore we can match each 
mode to the Bunch-Davies vacuum at a sufficiently large value of the physical frequency 
c s kHr e = M corresponding to the typical energy scale of the produced particles: 

%(r e ) ~ -r£f=ree- ic ^ , (3.18) 

which fixes 

The inhomogeneous contribution can be written as the integral of the Green's function 



fe £ (r) = -A-,/ VGM jWg-, , (3.20) 



with 

v-l 



G(t,t') = ^kc s r [Y u (kc s T)J u (kc s T')) - J v {kc s T)Y v (kc s r'\ . (3.21) 



In the limit r — > and r e — > — oo, the Green's function integral gives 



The final power spectrum therefore has two contributions, one proportional to 
the stochastic field and the other due to the quantum fluctuations of the field. 
These two contributions are uncorrelated, and so the power spectrum is the sum in 
quadrature of the two: 

_ (H\ * 2^I»^ , ( *m \ 2 <nH 

When the friction due to string production is a significant effect (A ^> 1), the 
first term in (3.23) (arising from de Sitter fluctuations) is very small, and the power 
spectrum is dominated by the second term. This was the regime considered in [27, 
28]. In a string theory realization, A is generically very small for unwinding inflation, 
but string production may or may not be the dominant source of perturbations 8z. 

As a check, (3.23) is derived by solving (3.13) directly for given 8n in Appendix 

D. 

Tensor power: The tensor power has a contribution from de Sitter perturbations 
(3.7). As mentioned above, string or particle production can also contribute to Vu [37], 
but we leave this for future work. 

Oscillations: Oscillations in the second term of (3.23) arise due to the periodicity 
of the sum (3.10). We estimate the amplitude of these oscillations relative to the 
time-averaged result (3.23) in Appendix E. We find that time averaging is a good 
approximation when HAt <C 1, where At = — because a high frequency driving 
force has very little effect on the mode. Instead, if HAt ~ 1, there are oscillations in 
the power spectrum due to string production with amplitude < 1. In the "prototype" 
model where At = l/2v ~ 1/2, the amplitude of the oscillations when HI — 1 is ~ 1CT 2 , 
and then falls off as a high power of HI. 

Tilt: When the de Sitter fluctuations in (3.23) dominate the power spectrum, the 
tilt is given by (3.9). If instead the fluctuations from string production dominate, the 
tilt arises from the second term in (3.23). To evaluate it requires knowledge of v, n, 
etc. In string theory these are calculable, and the results can be found in Appendix 
C. The average density of strings n ~ g(b,v)m z s / HI, where g(b,v) is a dimensionless 
function of the velocity v and impact parameter b, and (Hl)^ 1 is roughly the number of 
collisions per Hubble time. When b ^> m~ l , the stretched open strings are never light, 
and string production is exponentially suppressed. However, at least in the simplest 
cases one expects b to be small or zero. When b = 0, g ~ (9(10 _1 ) (Fig. 9). 

In general one expects a factor of H~ l in n, since the time-averaged density 
of produced particles is proportional to the number of collisions per Hubble time. 
This means that the second term in (3.23) varies during inflation only through its 
dependence on v. So long as the dependence of n on v is not too strong (as is the case 
in string theory for 0.5 < v < 0.999; see Fig. 9), its contribution to the tilt is small. 
The prefactor involving T(u) is very close to constant when A< 1. 
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To evaluate the tilt due to the factor of 7 3 , (3.11) implies 2:7 ~ 7 ~ V'/(6crH) 
when A < 1. At least in the "prototype" model H oc Q and V oc Q - 1/2 (2.7), 

7_ ^ Q 
7 Q(2Q-1)~2Q 2 " 

Since Q ~ Q — 2z/l, Q = —2z/l and hence 

dlwy -3 7 3i 1 

dink = ~ 3 H^/ = HIQ 2 ~ QiV' 

where we have used A ~ (Hl/2z)Q to derive the last expression. Since Q ^> 1 the tilt 
due to the factor 7 -3 is very small, and the spectrum of perturbations produced by 
strings is very close to scale invariant. 

This conclusion relied on a number of assumptions, and could be modified de- 
pending on the microscopic details of the model. For example, as mentioned above if 
a depends on the sizes of some compact dimensions and these change during inflation, 
there will be an additional contribution to the tilt. Another possibility is that the 
dependence of A,m , and n on v is not small and cannot be ignored. 

Still, at least in the examples we have investigated the tilt due to particle produc- 
tion is almost zero. This has the interesting consequence that there will be an "elbow" 
in the power spectrum: if the red-tilted de Sitter fluctuations dominate in the early 
phase of inflation when CMB perturbations freeze out (as is probably necessary for 
consistency with the observational constraints on n s ), then at some shorter scale the 
string fluctuations will surpass them in amplitude, and from then on the tilt will be 
very close to zero. Such a spectrum might have interesting observable consequences, 
for instance in so-called /i-distortions [39]. 

3.4 Non-Gaussianity 

The effective theory described by (3.2) will be non-Gaussian because of the DBI ki- 
netic term (the contribution from non-linearities in V is small [40]). The shape and 
amplitude of the non-Gaussianity from a DBI kinetic term was described in e.g. [41]. 
It is primarily equilateral: 

/inequilateral ~ (1 ~ <%) / <? 8 ™ . 

The observational constraints on DBI inflation from WMAP require c s > 0.054, or 
7 < 19. [41] 

If there are extra, approximately massless fields b, the situation is more complex. 
Assuming the fields b affect reheating - as they do in the case of unwinding inflation, 
since they parametrize the transverse separation of the brane and anti-brane that must 
annihilate to end inflation - they can reduce the amount of equilateral non-Gaussianity 
and increase the amount of local non-Gaussianity. This is described in for instance [42] 
and [43]. Roughly speaking, the bigger the relative contribution of 5b to (, the less 
equilateral non-Gaussianity there will be (because b are non-relativistic), and the more 
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local non-Gaussianity there will be. When (b) = at the time of reheating, the Z 2 
symmetry b — > — b prevents b fluctuations from affecting ( at least to lowest order. 
In this case neither the power spectrum not the bi-spectrum is affected by Sb, but 
higher correlations might be. On the other hand if (b) ^ 0, the Z 2 symmetry is broken 
and Sb is converted into ( with some efficiency, therefore increasing the level of local 
non-Gaussianity and decreasing equilateral. 

Another interesting feature that may arise in unwinding inflation is a peak at the 
"folded" shape. This occurs when there is particle production during inflation, or in 
general when inflation takes place away from the Bunch-Davies vacuum. If particle or 
string production is a significant contribution to folded non-Gaussianity could be 
a smoking gun indication of it [44]. 

Oscillations in the power spectrum of the type predicted by unwinding inflation 
can lead to unexpectedly large "resonant" non-Gaussianity [45]. Finally, if h(z) in 
(3.1) changes rapidly (as is the case for unwinding inflation on S 2 - see Appendix A), 
the non-adiabaticity may lead to additional non-Gaussianities. 

4 Stabilized compact ificat ions 

Unwinding inflation can occur in any model with at least one compact extra dimension 
threaded by a (p + 2) > 5-form flux F, and containing a brane charged under F. We 
have been assuming that the extra dimensions M. are stabilized by some effect other 
than the flux that is discharged during the cascade, and that the initial energy density 
in the unwinding flux is small enough relative to the stabilization mechanism that the 
geometry of M. is not strongly affected by the cascade. 

To show how this could happen, we give two examples. In both, gradually de- 
creasing the higher dimensional vacuum energy by discharging flux gradually decreases 
the 4D vacuum energy, as well as slightly decreasing the radius of the compact space. 
For a sufficiently slow decrease (combined with Hubble expansion to inflate away ki- 
netic energy), the radion field corresponding to the radius of the compact space should 
remain very close to its local minimum. 

(IS4 x Si: Stabilizing the Si can be accomplished with two ingredients: a positive 5D 
cosmo logical constant A 5 , and the Casimir energy of several bosons and fermions, at 
least some of which must be massive [34]. In our example, A 5 = A + fi 5 Q 2 /2 is a 
combination of a "bare" cosmological constant (which can arise from additional fluxes, 
vacuum energy, or other sources) and the flux that will discharge during the cascade. 

With these ingredients the effective potential for the radion of the Si can have a 
minimum, and the 4D vacuum energy A 4 (Q) in that minimum can be either positive, 
negative, or zero. If A 4 (Qo) > 0, the initial spacetime is metastable dS^ x Si. At some 
time, a bubble of brane appears and begins to discharge Q. For Q 3> 1, this is a 
gradual process that slowly reduces A 5 . Because the minimum of the radion potential 
occurs due to a balance of negative Casimir energy against positive vacuum energy, 
reducing Q has the effect of reducing A 4 and the radius of the circle. Numerically, in 
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the situation of interest the change in the radius of the circle is very small (Fig. 5) 
and can be neglected. 

Veff 




Figure 5. Stabilization of Si. Left panel: Effective potential for the radion field as the 
number of flux units decreases in the direction of the arrow. Right panel: Change in the 
minimum of the radion as the number of flux units varies. 

The weak point of this model is that Casimir energy is a quantum effect, while 
the energy in flux is classical. Hence this form of stabilization requires that quantum 
Casimir energy is larger than the energy in many units of classical flux. This objection 
is alleviated to some extent by the fact that the extra dimensions in our model need 
only be slightly larger than the string length, so that quantum effects can be significant. 
Furthermore in models with additional compact dimensions the flux quantization can 
be quite fine-grained, given only an 0(1) hierarchy in their volume relative to the 
string or Planck length. 

dSi x 5*2: The S 2 can be stabilized by positive 6D vacuum energy and magnetic 2- 
form flux F 2 that threads the S 2 [46, 47]. In this case the "negative" energy term in 
the effective potential for the radius of the S 2 is provided by its curvature. Once again, 
the 4D vacuum energy in the minimum can be positive, zero, or negative depending 
on the relative strengths of these contributions, and again, discharging part of the 6D 
vacuum energy by a flux cascade will gradually reduce the 4D vacuum energy. 

In this example, the flux cascade discharges a 6-form flux by the nucleation of a 
bubble of 4-brane, which is a circular string on the S 2 [3]. The string oscillates back 
and forth from pole to pole on the S 2 , wiping away an additional unit of 6-form flux 
with each pass. More details can be found in Appendix A. 

5 String theory 

The obvious context for unwinding inflation is string theory, which contains a variety of 
Z)-branes charged under Ramond-Ramond fluxes, as well as NS 5-branes charged under 
a 7-form flux. However, realizing any model of inflation in string theory is extremely 
challenging, because the compact dimensions must remain at least approximately stable 
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Figure 6. Stabilization of S2. Left panel: Effective potential for the radion field as the 
number of flux units decreases in the direction of the arrow. Right panel: Change in the 
minimum of the radion as the number of flux units varies. 

during 60 efolds of expansion of the non-compact dimensions. Furthermore, supersym- 
metry is broken at least at scale H, making the theory difficult to control. Here, we 
will only try to derive the scalings of the parameters of our 4D effective description 
with g s , m s , and the geometry of the extra dimensions. We will assume the compact 
geometry remains stable throughout the cascade, and ignore any couplings of the brane 
other than to the discharging flux. 

Unwinding inflation in string theory can in principle occur with a p-brane and 
flux F p+ 2 with 3 < p < 8. As we will see however, it appears to work most easily with 
p = 4 or 5, and is probably not possible when p — 8. 

5.1 Transverse volume 

In ten dimensions, there are 8 — p transverse dimensions that are not threaded by 
the unwinding flux. Classically, the brane bubble does not expand or extend in these 
directions, because there is no flux forcing it to. However, the relative position of the 
brane in the transverse directions (the impact parameter fields b) has a very strong 
effect on string production during brane collisions (see Appendix C), and in addition 
affects the time of reheating and potentially produces interesting non-Gaussianity (see 
Sec. 3.4). 

Furthermore, the size of the transverse dimensions determines the effective 4D 
brane charge. When they are larger than the string scale, the energy per unit flux 
is reduced by ratios of the string scale to the size of the transverse dimensions [6]. 
Because we require an initial flux number Qq > (9(100), the characteristic length of 
the transverse extra dimensions must be somewhat larger than the string length in 
order to avoid super-stringy or super-Planckian energy densities — although as we will 
see, a ratio ~ 5 suffices for p = 5. The lack of transverse dimensions appears to rule 
out p = 8, and makes p = 7 problematic. Perhaps this conclusion can be modified by 
considering warping, but we will not do so here. 
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5.2 The flux cycle 

The flux F threads a non-contractible (p — 2) cycle in Ai. If the initial radius Rq of 
the brane bubble when it nucleates is small compared to the size of that cycle, the 
bubble will be roughly spherical and centered at a point on the cycle the flux threads. 
After nucleating it will expand in all the flux directions and collide with itself as it 
wraps around, initiating a flux cascade. The geometry of the cycle is important in a 
number of ways, affecting both the average rate of discharge of background flux and 
the details of the perturbations (Appendix A, Appendix B). There are too many such 
possibilities to discuss exhaustively here, so we will focus on some specific cases (see 
[3] for some additional examples). 

If the would-be radius of the bubble Ro is larger than all the length-scales of the 
cycle, the dominant instanton will not have spherical topology or be localized at a 
point in the cycle. Instead, it will be a sphere in the 3+1 non-compact dimensions of 
spacetime, but will wrap the cycle. In this case there is no flux cascade — the bubble 
can be described by dimensional reduction, and so is just an ordinary 3+1D bubble 
inside of which the flux is reduced by only one unit. 

If the cycle is anisotropic, the brane bubble may wrap a smaller sub-cycle but 
be localized and expand in other, larger directions [23, 48]. For instance, consider a 
cycle of the form Si x C. If I is the radius of the Si and d is the characteristic length 
scale of C, when I > Rq > d the dominant instanton should wrap C and have spherical 
topology in the dS± x Si directions. In this case the dynamics can be described by 
dimensional reduction on C, resulting in the "prototype" model of Sec. 2.1 with an 
effective 3-brane with tension a ~ a p d p ~ 3 . 

Even if the cycle is not anisotropic and the initial bubble is small, it is possible 
that brane interactions could lead to an attractor solution where the brane wraps some 
subcycles. For example, as the brane freely expands as a sphere on a torus T q it winds 
around more and more densely (its overall length increases as z q ~ l ). Eventually, even 
at high velocity brane interactions may allow reconnection to a configuration with few 
or no self-intersections — namely a configuration where the brane wraps all but a single 
Si of the T q . If this quasi-equilibrium is reached before 60 efolds from the end of 
inflation, the model will be effectively described as above. A few details of unwinding 
inflation on T q may be found in Appendix B. 

Another interesting case is that of p = 4 where the 6- form flux extends in dS± x S^. 
As described above, at least from the 6D point of view the S2 can be stabilized by 2- 
form flux and a 6D vacuum energy, of which part is the 6-form flux F that discharges 
during unwinding inflation. The bubble of 4-brane is an oriented string on the S2 
that appears as a circle around some point. The circular string expands, reaches its 
maximum extent at the equator of the S2, and then contracts on the antipodal point 
from the one where it nucleated. Assuming it does not annihilate, it will invert, and 
then again expand until it contracts on the point where it initially appeared. Each 
such cycle discharges one additional unit of flux. The 4D effective action for the brane 
and its perturbations is described in Appendix A. 
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5.3 Effective parameters from string theory 

In this subsection we will relate the parameters of our 4D description to those of 
the underlying string theory. We will work in conventions where m 2 s = l/(2ira') = 
l/(27rZg). The fundamental tension and charge of a p-brane of type IIA/B string 
theory are 

m P+i 

a n - 



p g s {2Tr)(p-W 

When the p-brane wraps a p — 3 cycle with volume d p ~ 3 , the tension of the resulting 
effective 3-brane is 

m p+l d p-3 



a = a n d p 



-3 . "°s 



p g s {2Tr)(p-W 

When all compact dimensions have radius ~ d, except an Si which has radius 
the 5D effective charge after reducing on everything but the Si is 

5 aM^~ M d 2 ^ _ (dm s ) 2 Pg^ p - 3)/2 
^ d 5 ~ (2vr) 2 rf 11 m6 ' 

The Hubble constant during inflation is 

tt2 87tG 4 /i 5 Q 2 / 

H ~ . 

3 2 

In order for the flux cascade to occur, we need R = 4cr / 'n < 1/2, so that the 
dominant instanton is localized on the Si (as opposed to wrapping it). For simplicity 
we will require Rq > d so that the brane wraps all but the Si of length I. Furthermore 
we require that Hd, HI < 1, on the prejudice that de Sitter compactifications with 
d, I > 1/H are difficult or perhaps impossible to stabilize. We require /, d to be larger 
than the string length, and N > 60 efolds of slow-roll inflation to solve the curvature 
problem. For consistency with observational constraints we require that the scalar 
density perturbations have the observationally correct magnitude. 

We can satisfy all these constraints simultaneously for 4 < p < 6. Qualitatively, 
p = 8 is not possible because there are no transverse dimensions to dilute the energy 
per unit flux, and so one cannot have Qq ^> 1. Similarly, the case p = 7 is incompatible 
with these constraints because the single transverse dimension with length d should 
satisfy Hd < 1, while simultaneously being large enough to dilute the flux. 

For p = 6 the Lorentz factor is typically large (7 ~ 50), and the parameters would 
require some tuning for consistency with the current constraints on non-gaussianity 
[41]. 

For p — 3 all the constraints above can be satisfied, but at the cost of a large 
hierarchy l/d ~ 40, and string production — which has a nearly flat spectrum — is the 
dominant source of perturbations. 

Examples: The two examples that work the most easily are p = 4 and p = 5. 
We give a table (Table 1) with a set of parameters for each case, and the resulting 
observational parameters. In both cases the scalar spectrum is dominated by de Sitter 
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p 



9, = 0.01 


V c = 2.4 x 10~ 9 


I = 20m; 1 


V h = 5.0 x HT 11 


d = 2.0m- 1 


r = 2.1 x 10~ 2 




n s -l = -0.032 


Qo = 400 


H = 0.05m s 


g* = 304 


7 = 11.5 



g s = 0.05 


P c = 2.4 x 10~ 9 


/ = 20m7 1 


V h = 1.4 x HT 11 


d = 4.9m;: 1 


r = 5.8 x lO- 2 


6^0 


n s - 1 = -0.032 


Qo = 400 


# = 0.04m s 


g* = 314 


7 = 22.9 



Table 1. Parameter sets and the corresponding cosmological observables for p = 4 and p = 5 
branes that expand around an S\ with circumference I, and wrap one or two directions of a 
torus with circumference d. Both sets satisfy all the constraints mentioned in the text. 



fluctuations rather than string production for 6^0: strings account respectively for 
10% of the power in the p = 4 case and 1% for p — 5. Increasing the impact parameter 
b suppresses the string contribution even more. 

In both of the above cases, the oscillations in the energy density V{z) are (D(Q~ 2 ) ~ 
10~ 5 , too small to be observed in the power spectrum. However the oscillations in the 
perturbations due to string production are relatively large (see Appendix E). If strings 
account for 10% of the scalar power as in the example above with p = 4, this could 
produce observable oscillations in the power spectrum. Another interesting feature is 
the presence of an "elbow" in the power spectrum, where the perturbations due to 
strings overtake the red-tilted de Sitter perturbations at high k. 

6 Reheating 

When a brane and anti-brane pass within a string length of one another, the lowest lying 
mode of the open strings stretched between them is tachyonic and will spontaneously 
condense. If the relative velocity of the brane /anti-brane pair is relativistic, the branes 
will pass by each other and separate to more than a string length, in which case the 
lowest stretched string mode is no longer tachyonic (and one can meaningfully speak 
of its number density). The calculation in Appendix C shows that the density of such 
"tachyons" is small in string units when the velocity v > .5. The conclusion is that 
the tachyon does not have time to condense when the velocity is relativistic; the same 
conclusion was reached by different methods in [3]. 

If instead the relative motion is slow and the brane and anti-brane spend a sig- 
nificant amount of time separated by less than a string length, the tachyon can fully 
condense. This is interpreted as the annihilation of the brane and anti-brane, with the 
~ la worth of energy density converted into closed string radiation. Once the tachyon 
has condensed the open string physics is strongly coupled [49], and it is believed that 
the would-be open strings ending on the branes get confined into closed strings. [50] 

Other models of brane inflation use tachyon condensation as a mechanism for 
reheating, but unwinding inflation is unusual in that the brane /anti-brane pass close 
together many times before the tachyon condenses. At first, the flux is large enough to 
keep the branes moving relativistically. Eventually, however, the flux decreases to the 
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point that the branes move with non-relativistic velocity. At some point after this, it 
will annihilate with an image anti-brane, reheating the universe and ending inflation. 
This may occur at Q — 0, or at some Q e satisfying Q ^> Q e > 0. 

A potential problem that could arise is if the brane sometimes over- or under- 
shoots Q e in some Hubble volumes, meaning that there are regions where the number 
of flux units remaining differs from Q e . These regions will be surrounded by tensionful 
brane. If they are sufficiently small, the tension of the brane will cause them to 
collapse and annihilate into radiation (or perhaps form black holes that subsequently 
evaporate). But if they are sufficiently large and the energy density in them is lower 
than their surroundings, they may (depending on the tension, see e.g. [51]) expand. If 
so, they will collide and percolate, replacing the surrounding phase with the one inside 
them. 

This process will create a feature in the power spectrum at the length scale cor- 
responding to the typical separation between such "overshoot" regions. Unless this 
length scale is ~ e N */H e (where H e is the Hubble scale at the end of inflation) it is 
unlikely to affect CMB perturbations, but might have other interesting consequences. 
Like percolating bubble collisions in first-order phase transitions, it could create grav- 
ity waves [52] or primordial black holes. A detailed analysis of the model is required 
to determine this. 

If the 4D vacuum energy corresponding to Q = is negative, then overshoot 
regions with sufficiently small or zero Q may have negative vacuum energy. If these 
regions expand, they may crunch the universe. Presumably, we live in the vacuum 
with the lowest Q accessed by unwinding inflation. 

Models with very high re-heat temperature may suffer from a monopole problem, 
although the details of this are strongly model dependent. Reheating may also produce 
cosmic strings, with potentially observable consequences. 

7 Conclusions 

From the "bottom-up" or effective field theory point of view, inflation should be de- 
scribed by the simplest model consistent with data — perhaps single-field m 2 2 . On the 
other hand, natural phenomena sometimes turn out to be described by field theories 
that are much more complex than initially seemed necessary, with the Standard Model 
perhaps the prime example. From the low-energy effective point of view, unwinding 
inflation looks like an unholy union of m 2 </> 2 , DBI inflation, trapped/dissipative infla- 
tion, hybrid inflation, and contains multiple light fields, an oscillating potential that 
may lead to resonant non-Gaussianity, and particle and graviton production. But from 
the microscopic point of view it arises from a very simple mechanism and seems to fit 
naturally into string theory, itself a theory based on one very simple assumption. 

Many issues remain largely unexplored. The dynamics of re-heating should be 
more carefully investigated. The stability of the model — both of the compact extra 
dimensions and to other potential instabilities — needs to be studied. If the formation 
of the initial bubble is not too rare, collisions between two unwinding bubbles are 
possible, and their dynamics would be interesting to investigate. The effects of the 
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geometry of the compactification manifold on the power spectrum of perturbations have 
proven rich and interesting already at the superficial level we have investigated them. 
One particularly important issue is the effect of the light scalars corresponding to the 
transverse directions on the fluctuations (both their amplitude and non-Gaussianity). 
The spectrum of tensor modes and amount of dissipation from brane Bremsstrahlung 
remains to be calculated. Realizing unwinding inflation in string theory is a very 
important goal in its own right, and would help frame all of these questions. 
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A Unwinding inflation on a sphere 

In this section we consider a 4-brane on dS± x S2 with an initial six-form field strength 
F e = fi 6 / 2 Qo (where fj, is a parameter with dimensions of mass, and Q counts the units 
of flux). The radius R of the S 2 can be stabilized by 2- form flux; see Sec. 4 for details. 
We use the coordinates: 

ds 2 = -dt 2 + e 2Ht d* 2 + R 2 (d6 2 + sin 2 9d<j) 2 ) (A.l) 

and ignore the negative spatial curvature in the de Sitter directions. 

The cascade will begin with the nucleation of a 4-brane extended in the 3 spatial 
de Sitter dimensions and co-dimension 1 on the sphere. On the sphere, the brane will 
be an oriented, circular string with angular radius 6b(t = 0) centered on a random point 
that we choose to be 9 = 0. When the bubble first forms it will contain Q — 1 units of 
flux, but it will expand until it reaches the equator, contract on the point 9 = tt, and 
then — assuming it doesn't annihilate — invert and expand, with opposite orientation, 
back towards the equator, contract to 9 = 0, etc. It will continue to oscillate between 
9 = and 9 = n, discharging one unit of flux on each pass, until the cascade ends by 
annihilation. 

The brane/flux action for this configuration is: 

S = J d 4 xd9d(pe 3Ht R 2 sin 9 x 

W - - R 2 (V - e-^(ve b y - j^mr) - ■ 

(A.2) 
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Background To solve for the background evolution, we treat the position of the 
brane as a function of time only: 9b = 8b(t). The action is: 



S 



d x 



e 3Ht R 2 sin 9 [ -aS(R(9 - 9 h ))\l I - R 2 9, 



• 2 



2 • 6! 



(A.3) 



It is useful to introduce a monotonically increasing variable Ub to keep track of 
the number of passes that the brane has made over the sphere: u>b = 9b for the first 
pass, Ub = 2n — 9 b on the second, u b = 2n + 9 b on the third, and so on. Using this 
variable, the energy density in the flux is: 



F 2 M 



Q 







LOb 
7T J 



Qo 



LOb 
L 71 



9 ((-1) (COS £J 6 - COsfl)) + 

cosub + cos 9)) 



(A.4) 



where, again, [. . .] denotes the integer part. 

Performing the integrations over and 9 and defining z = Rub gives 



S 4 



d A x2rre zm 



-Ra\ sin — |vT 
R 



V A {z) 



(A.5) 



with 



jx e R 2 
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^R 2 



Qo 
Qo 



Z ' 
7tR. 

z ■ 

7rR. 



1M i 



■1) 



[z/{<kR)\ 



COS 



R 



Wi + (_1)[*/(^)1 



cos 



R 



(A.6) 



(see Fig. 7). The background equation of motion is 



z + 



3H 



-z + 



cot 



+ 



i?7 2 i?0"7 3 1 sin 



(A.7) 



R\ 



where 7 = (1 - i 2 )" 1 / 2 and V{{z) = —fJ, 6 R\ sin ^|(Qq - 1/2- [^]) = -ft|sin||. 

The force term proportional to V is the electric force, which pushes z to increase 
until the flux is discharged and Qo — 1/2 — [-^] reaches zero or changes sign. The 
curvature of the sphere induces an additional force term cot (z / R) / (R-y 2 ) that changes 
sign at the equator and has a possible divergence at the poles. This arises from the 
fact that the brane's area is maximized at the equator and zero at the two poles. 

We have not been able to solve (A.7) analytically, although it can be solved in 
terms of Jacobi functions when H — V — 0. These solutions (which describe a free 
string oscillating on the sphere without Hubble friction) are a good approximation near 
the inversion points z = nnR, and demonstrate that the kinetic energy 7I sin(z/R)\ is 
finite and non-zero there, indicating that 7 ~ 1/| sin(z/R)\ near the inversion points. 
Equation (A.7) can be solved numerically without difficulty (Fig. 8). 
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Perturbations To determine the perturbations around the smooth background evo- 
lution we define 9 b (t,x,(f>) = 9 h {t) + 59 b [t ) x ) (f)) and expand (A. 2) to second order in 
56b- After performing the 9 integration, the action for the perturbation is: 



S 



d xa<pe 
1 



3Ht 



-R 3 a sin 9 b ^{59 b f + R 3 a cos 9 b j9 b 89 b 59 b 



— -R ersin^Vy 



(d 3 69 b ) 

c 2Ht 



-Ra^d^bf + 



Ro_ sin 9 b _ 
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Expanding the dependence of the perturbation in modes: 



5z = J2 Sz n e in<t> , 



591 



(At 



where z is defined as above. Performing the <f> integration gives: 



d 4 xire 3Ht 



sin 



R 



Ra5z n + 2<T cot — z5z n 5z n — Ra^ ^jr — 

R e zttt 



+ 



a 

fry 



a'yn 



V"{z) 



R sin 



5zl 



sin 



R\ 



Dropping the subscript, the equation of motion is: 

'cot 



5z + 



3H + z 
-1 

R 2 sin 2 | 7 2 



i? 



5z + 



+ 



A; 2 



+ 



7 2 e 2ffi ' i? 2 7 2 Isin 



+ 



Hi 



i?u7 3 dz \ sin 



= 0. 



(A.9) 



(A.10) 



This equation is somewhat difficult to deal with numerically, primarily due to the co- 
efficient of 5z. However, at superhorizon wavelengths (k = 0), 5z = z is an exact 




Figure 7. Left panel: The pressure V for the compactification on S%. Right panel: The 
potential V around the minimum; the smooth line is the quadratic approximation. 
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Figure 8. From left to right: the "unwound" angular position of the brane oj^ = z/R, its 
angular velocity z/R, and its Lorentz factor 7 multiplied by | s'm(z/R)\. 

solution to (A. 10). Because the time-averaged friction (the time average of the coef- 
ficient of 5z) is positive, this should be the "growing" solution, a conclusion we have 
verified numerically. Since z is a constant plus small oscillations, the perturbations are 
nearly constant outside the horizon. To determine their amplitude requires numeri- 
cally solving (A. 10) for finite k. Preliminary investigations indicate that their average 
amplitude is parametrically similar to the S\ case (3.6), but with additional oscillations 
with amplitude depending on HR. 

Calculating the degree of string production when the brane turns inside out is 
an interesting problem that we leave for the future. We tentatively expect that the 
results will again not be substantially different from that of parallel brane/anti-brane 
scattering, but this remains to be established. 

B Unwinding inflation on a torus 

Here we consider the case of a flat g-torus, where the bubble nucleates with a diameter 
smaller than any cycle of the torus. Ignoring interactions for now, the brane will expand 
freely and uniformly in all q directions with velocity v ~ 1, wrapping repeatedly around 
the cycles of the torus and intersecting itself. Its overall area will increase as f2 9 _iz 9-1 , 
the area of a q — 1 sphere of radius z. 

The flux will discharge in a complicated pattern of overlapping spheres, but the 
average number of flux units discharged is simple to calculate. In the covering lattice 
of the torus, when the radius is z the wall of any image bubble that nucleated within 
a distance z of the origin will have crossed the origin and discharged a unit of flux. 
Since there are ~ (Q q ^i/q)(z/l) q such points, the flux 

Q(z,y) = Q -^(^)\l + r(z,y)), 

where y are the coordinates on the T q , l q is the volume of the torus, and r(z, y) is 
the error in approximating the number of points within a sphere by the volume of 
the sphere over the volume of the torus. The function r(l,y) ~ 0(1), it oscillates 
and falls off like a negative power of z/l for z ^> I and it averages to zero over y : 
J d q yr(z,y) = 0. 
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Dropping the small oscillating remainder term r, the action is 



J dtdxe 3Ht i^ q+2 n q ^- l Vi^¥ + (q - ^ (y)") 2 } ' (R1) 



The equation of motion is 



3Hz q-1 V 



7 2 7^ ^Vt q _iaz q - 
The power spectrum of de Sitter fluctuations is (3.6) 

# 4 



c 



ATi 2 VL q _ l a q+2 v 2 z<i- 1 ' 



Physically, the az q ~ l in the denominator represents the fact that the effective mass of 
the brane is increasing as it expands on the torus. The larger the mass, the greater 
the inertia and the smaller the amplitude of de Sitter fluctuations. 
The tilt is 

(q-l)v A H 

n — 1 w — — 1- 4 — 

s Hz + # 2 ' 

which is red and relatively large in magnitude. 



C String production from brane/anti-brane scattering 

In this appendix we derive the amount of open string production between a parallel 
p-brane and anti-p-brane that pass each other with impact parameter b and velocity 
v . 5 After a few efolds of unwinding inflation on Si the radius of the brane bubble 
is exponentially large, so that parallel branes are a good approximation (see Fig. 2). 
We perform the calculation in flat space, ignoring the de Sitter curvature — but since 
string production takes place in a period l s /v <C H~ l , we expect this is also a good 
approximation. 

We ignore closed string production, which may occur via brane Bremsstrahlung or 
by decay of the open strings that are produced. At ultrarelativistic velocities and large 
brane co-dimension, closed string radiation is a very important effect 6 and neglecting 
it may not be justified. However, Bremsstrahlung is proportional to acceleration, and 
therefore should not affect the motion of the brane at its terminal velocity (2.4). In 
principle stretched branes could also be produced, but the rate should be suppressed 
relative to open string production 

At lowest order the rate of string production in brane scattering is determined by 
the imaginary part of the annulus diagram with the appropriate boundary conditions 
[36]. Physically, string production occurs because the mass of the stretched strings 
has a contribution that is proportional to length, which is changing with time as the 

5 A calculation using effective field theory in a related scenario appears in [53]. 
6 We thank L. McAllister for discussions on this point. 
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branes move. The calculation is T-dual to the pair production of open strings in an 
electric field, [54] and has many features in common with the Schwinger effect. 

Our analysis differs from [36] in that we consider brane/anti-brane scattering, 
rather than brane/brane. This has important consequences at low velocity: the 
brane/brane case becomes supersymmetric as v — > 0, and the rate of string production 
scales as v to a positive power. This is related to the fact that the lightest stretched 
open string mode is massless when the branes coincide. For brane/anti-brane scatter- 
ing, there is a tachyonic mode that will condense if the brane and anti-brane remain in 
close proximity for long enough. As we will see, relativistic velocity prevents this from 
happening even when the impact parameter 6 = 0: the lightest mode is tachyonic only 
briefly, and does not condense. Instead, a finite density of strings are created in this 
mode. 

Our starting point is equation (11) of [36] for the phase shift for forward scattering: 
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(CI) 

where V p is the volume of the p-brane (or anti-brane), #'s and r] are the Jacobi and 
Dedekind functions for the definitions and properties of which we address the reader 
to [55], x — tanh _1 t> = tanh _1 t> 2 — tanh -1 v\ is the relative rapidity, and b = \b\ is 
the impact parameter. The charges (e a ) in the brane-anti brane case are such that 
e 2 = -e 3 = -e 4 = -1. 

We are interested in the imaginary part of (C.l). The poles of the integrand 
are the zeros of 9i(xt/2ir\it/2), i.e. x^l^ = k G Z . The residue of the k-th pole is 
l) k . Therefore: 
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which, using properties of the Jacobi functions, can be rewritten as 



X 
(C.2) 
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(C.3) 



We want to expand the term in the curly brackets around small q = exp(2iri(i7rk/x)) 
to reveal the string mode expansion. One has 



V- 12 (q) [("I + ("1)% 4 (?) + (1 + (-1)^4 (?)] 



-16+72(-l) fc v /^-256g+0(g 3/2 ). 

(C.4) 



The integer coefficients are the number of stretched open string degrees of freedom 
at the corresponding oscillator level — a complex scalar tachyon, 16 massless fermionic 
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degrees of freedom, etc. — while the factors of (— l) fc arise due to spin statistics as usual 
for the Schwinger effect (c.f. [54]). 

Ordinarily the sum over k is almost irrelevant, because the higher terms are 
exponentially suppressed. In our case the sum looks problematic due to the term that 
corresponds to the tachyon (at least for small h where the string is actually tachyonic 
at closest approach): 



^(-l^AT^^exp 



k=l 



X 



However, even for 6 = the sum can be performed exactly in terms of a polylogarithm, 
and (with the help of the (— l) fc factor) is real and finite away from v = x — 0. The 
physical interpretation is that when v ~ 1 the stretched string is only very briefly 
tachyonic and does not have time to condense. By contrast for v <C 1 there is a 
divergence, which leads to brane-antibrane annihilation due to tachyon condensation. 

From (C.3), the amplitude given by the four lightest open string modes (tachyonic, 
massless and the first two massive) is: 



T lm p s V p , 2 
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The energy density in the produced strings is simply related to (C.5): we should 
multiply by the energy per string and divide by the volume. The energy per string is 
evidently [56] 



E 2 = v - [mt(z 2 + b 2 ) + 7rm 2 s (n - 1)1 « - 2 m 4 s (z 2 + b 2 ), (C.6) 
X X 

where \/ z 2 + b 2 is the length of the string and the integer n is the oscillator number 
(in our notation n = is the lowest mode of the stretched string, which is tachyonic 
when z = b = 0). The pre-factor v/x can be derived by matching to a field theory 
computation [56]. In any case (because x ~ l n 7) it is not an important effect except 
for ultra-relativistic scattering, and reduces to the standard form when v <C 1. The 
approximation in the last step holds for unwinding inflation, where the typical value 
of z ^> m" 1 and only the first few oscillator modes are produced. 

The energy density of strings produced at a single collision is then 

Pc ( x ,z) = EIm(6)/V p *±^ p x p/2 ^F(b,x)^+¥. (C.7) 

For 0.9 < v < 0.999 and 6 ^ the quantity \ j^yX p/2 ^F ~ (2vr) 3 " p , as is its derivative 
with respect to v (Fig. 9). For g s 1, this means that the energy density in produced 
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Figure 9. The quantity n c m^j (2ir 3 ~ p ) (solid line) and its derivative with respect to v (dashed 
line), in units m s = (2irl 2 )~ l > 2 = 1. Both are 0(1), with only a mild dependence on v for 
0.5 < v < 0.999. The case shown is p = 4, but the magnitude and shape are very similar for 
other p. The divergence as v — > is due to the tachyon; the divergence as v — > 1 is due to 
the copious production of strings at ultrarelativistic velocities. 



strings is much less than the brane tension a ~ m p s +1 / g s at least for z ~ l s . For 
z 3> l s the energy density in stretched strings will exceed the brane tension, but 
during unwinding inflation the brane continues to move at nearly constant v due to 
Hubble dilution of the strings and the force from the background flux; see Sec. 3.3. 
In the notation used in Sec. 3.3: 

YJlP 

2H ^ 3 1 3HAt V dz 7 3HI V dz 1 

D Effect of string/particle production on perturbations 

In Sec. 3.3, we solved (3.13) using a Greens' function method. Here, we will solve it by 
another method that serves as a check and adds some intuition regarding the behavior 
of the solutions. From (3.13) and (3.16), the equation we wish to solve is 




Recalling that an overbar denotes time-average, and that df /dz is constant, this equa- 
tion describes a damped harmonic oscillator with a spring constant that decreases with 
time, and a constant driving force that also decreases — but more slowly than the spring 
constant. The oscillator is underdamped at early times, but becomes overdamped at 
a certain critical time. 
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Switching to conformal time, we have 
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where A = (df/dz)/2Ha'y 3 , and B = -Wq \J Hn / 2o"7 3 are dimensionless parameters 
related to the amount of string or particle production. 
The general solution is 
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where ft is the regularized hypergeometric function. 

The late-time (r — > 0) limit of this expression is simple: 



2 2+v T(l + v) ny 



5z k -> -C 2 



To fix the coefficient C*2, we expand (D.3) in the limit r — > — oo, and match it to the 
Bunch-Davies state 

ziZV 



-ifcr/7 



2y^fk 

(see (3.18)) at r e = jM/k. Solving gives 
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2H f H\ 



A/2 



2 -A/2 



r(l) 



where 0i 5 2 are real. Because the phase of B is random and uncorrelated with the 
first term, the power spectrum will be the sum in quadrature of these two terms, and 
therefore 4>i,2 drop out of the result. The late-time perturbation is 



Szk ->■ 



V2T{u+l) 

k 3/2 

in agreement with (3.23). 
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dz dz dz 




Figure 10. Left panel: the behavior of the modes bz in the presence of particle production, 
when the friction term A > 1. The additional friction damps the oscillations rapidly, and 
after they become overdamped the source term B pushes bz to a non-zero value. Center 
panel: A < 1 and B < H 2 j^fa. Right panel: A < 1 and B > H 2 j^fo. 



The solutions (D.3)have three regimes of interest (see Fig. 10). When A ^> 1 the 
mode is more highly damped than is the case from Hubble friction alone, and therefore 
the amplitude of the homogeneous modes is exponentially supressed. However, the 
source term B may compensate. When A < 1, the damping from string/particle 
production is negligible compared to Hubble friction, and the homogeneous modes are 
undamped. However, there are still two regimes: B ^> H 2 /y/cr, in which case the 
fluctuations due to string/particle production are larger than those produced by de 
Sitter fluctuations, or B <C H 2 in which case the contribution of string/particle 
production is negligible. 



E The amplitude of oscillations from string or particle pro- 
duction 

Here, we estimate the error made in time averaging the equation for the perturbations 
due to string or particle production. The equation for the perturbation is 

2^a8z + QHadz-f + V"{z)8z + 6f = 0, 

where 

i i v v ' 

din din 
5U-gf n c (3H Q(t - U) - 8(t - U)) + Sn c -^ B(t - t t 

When the number of collisions per Hubble time is large, one expects time averaging 
to be a good approximation. But we may be interested in scenarios where the number 
of collisions per Hubble time is order 1. 

The terms proportional to 5z are friction terms for the perturbation, but their 
average value will turn out to be much less than H (at least in string theory), and 
therefore they have little effect with or without averaging. The term 5ti(3H Q(t — tj) — 
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5(t — ti)) oscillates and cancels after time averaging. We can estimate the amplitude 
of the oscillations caused by these terms as follows. The effect of the S function is 
cancelled by the function on average. Therefore the amplitude of the oscillation is 
the "kick" provided by the 5 function. The equation for the perturbation near the time 
of a collision will be approximately 

o... dm c dm c . 

2<T7 dz = -q— n c d(t - t^bU = n c b(t - U)bz/z. 

Integrating both sides gives 

CfTH 

Abz/bz = ^n c /(z2a 1 3 ), 

and therfore 

(Iff) 

Abz/bz = Abz/Hbz = -^n c /{z2Ha^). 

This can be evaluated exactly given the underlying model. To get a rough sense of its 
magnitude, if n c ~ mf, dm c /dz ~ m 2 s) and a ~ m\/g s , one obtains 

Adz /dz ~ „ . 

This is very small in the models considered in Sec. 5. 

Lastly, the term proportional to 5n c is a driving force that displaces the mode from 
its minimum at 5z = 0. The combination of the G functions and the decaying expo- 
nentials make the sum oscillate around its average value, which is ~ 5n c ^ £ /(3iJAt), 
where At = ti — tj_i. The amplitude of the oscillations in the driving force is simply 
bnj^f, and therefore the amplitude divided by the average is 

A5n c /dn c ~ 3HAt. 

The perturbation 5z due to string/particle production is proportional to the driving 
force (see (D.4)) when the driving force is constant. Therefore oscillations in the 
driving force with frequency u> = 2ti / At < H should introduce oscillations in 5z with 
the same amplitude. Hence when HAt ~ 2ir, the oscillations in the power spectrum 
due to string production are ~ 1. 

We have verified this conclusion numerically, and found that the amplitude of the 
oscillations is 10~ 2 when HAt ~ .5, and fall off thereafter like a high power of HAt. In 
the models of Sec. 5, this is the leading source of oscillations in the power spectrum. 
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